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ABSTRACT 

We study the X-ray emission of baryon fluid in the universe using the 
WIGEON cosmological hydrodynamic simulations. It has been revealed that 
cosmic baryon fluid in the nonlinear regime behaves like Burgers turbulence, i.e. 
the fluid field consists of shocks. Like turbulence in incompressible fluid, the 
Burgers turbulence plays an important role in converting the kinetic energy of 
the fluid to thermal energy and heats the gas. We show that the simulation 
sample of the ACDM model without adding extra heating sources can fit well 
the observed distributions of X-ray luminosity versus temperature (L x vs. T) of 
galaxy groups and is also consistent with the distributions of X-ray luminosity 
versus velocity dispersion (L x vs. a). Because the baryonic gas is multiphase, 
the L x — T and L x — a distributions are significantly scattered. If we describe 
the relationships by power laws L x oc T° LT and L x oc a aLV , we find > 2.5 
and o>lv > 2.1. The X-ray background in the soft 0.5 — 2 keV band emitted 
by the baryonic gas in the temperature range 10 5 < T < 10 7 K has also been 
calculated. We show that of the total background, (1) no more than 2% comes 
from the region with temperature less than 10 6 5 K, and (2) no more than 7% is 
from the region of dark matter with mass density pdm < 50pdm- The region of 
Pdm > 50pdm is generally clustered and discretely distributed. Therefore, almost 
all of the soft X-ray background comes from clustered sources, and the contribu- 
tion from truly diffuse gas is probably negligible. This point agrees with current 
X-ray observations. 
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1. Introduction 

In the linear regime of gravitational clustering of cosmic matter, the evolution of baryonic 
gas follows dark matter point by point (Bi et al. 1992; Fang et al. 1993; Nusser & Haehnelt 
1999; Nusser 2000). That is, the density and velocity distributions of baryonic matter can be 
obtained by a similar mapping from the dark matter field. The similar mapping has also been 
used in modeling the gas in clusters of galaxies (Kaiser 1986). However, it has been found 
that the similarity between the mass and velocity fields of baryonic gas and dark matter 
is broken in clusters that are formed through highly non-linear evolution. The mechanism 
of the breaking of similarity is due to the different dynamical behavior of dark matter and 
baryonic gas. The former is collisionless, while the latter is approximately an ideal fluid. 
The velocity field of collisionless dark matter particles is multivalued at the intersection of 
their trajectories, while baryon fluid always has a single-value velocity field. Thus, shocks in 
baryon fluid will appear at the intersection and break the similarity between the mass and 
velocity fields of baryon fluid and dark matter (Shandarin & Zel'dovich 1989). 

At later times, it has been recognized that shocks occur not only in high-, but also in 
middle- and even low-density regions. This point has been shown by the dynamical equa- 
tion of baryonic gas. Although cosmic baryonic gas is a Navier-Stokes fluid, its evolution is 
dominated by the growth mode, which is approximately governed by a random-force-driven 
Burgers equation, and the random force is produced by the gravity of the random field 
of dark matter (Gurbatov et al 1989; Berera & Fang 1994; Vergassola et al. 1994; Jones 
1999; Matarrese and Mohayaee 2002; Pando et al. 2002, Pando et al 2004). When the 
Reynolds number is large, shock-caused turbulence, called Burgers turbulence, will develop 
in the fluid (Lassig 2000). The dynamics and thermodynamics of baryonic gas will be sub- 
stantially affected by the Burgers turbulence. Dynamically, it will lead to the discrepancy 
between baryonic matter and dark matter, like the discrepancy of a passive substance from 
the underlying field during nonlinear evolution (e.g. Shraiman & Siggia 2000). Thermody- 
namically, the Burgers turbulence leads to the multiple phases of thermal properties and to 
converting the bulk kinetic energy of the fluid to thermal energy. For the system of cos- 
mic baryonic and dark matter, a large Reynolds number actually corresponds to the onset 
of the non-linear regime of the gravitational clustering. Therefore, the similarity between 
baryonic and dark matter will inevitably break in the non-linear regime, and the dynamical 



- 3- 



and thermodynamical effects of the Burgers turbulence must be considered. 

Some effects of the Burgers turbulence have been detected with observational data 
and/or hydrodynamic simulation samples of cosmic baryonic gas. First, the statistical de- 
coupling between the density and velocity fields of baryonic gas and dark matter is found 
to be significant on scales larger than the Jeans length (Pando et al 2004; Kim et al 2005). 
The importance of the discrepancy in the early nonlinear evolution has also been noted by 
Yoshida et al. (2003). Second, the transmitted flux of QSOs' Lja absorption is found to 
be remarkably intermittent (Jamkhedkar et al 2000, 2003, 2005; Pando et al 2002; Feng 
et al 2003), which is consistent with the prediction of the intermittence of fully developed 
Burgers turbulence (Polyakov 1995; Balkovsky et al. 1997; Frisch et al 2001). Third, the 
temperature field of the baryonic gas is multiphase, and the heating caused by the shocks of 
the Burgers turbulence is substantial. Consequently, the distribution of the baryon fraction 
on large scales is nonuniform (He et al. 2005), and high-entropy gas is produced (He et al. 
2004). 

In this paper, we extend these studies to the X-ray emission of the hot baryonic gas 
in the universe. X-ray emission of galaxy groups is considered as important evidence of 
the similarity breaking between baryon fluid and dark matter. The relation between X-ray 
luminosity L x and temperature T predicted by the similarity takes the form of a power law, 
L x oc T a , with an index of a = 2 (Kaiser 1991), while the observed result is a > 2 (Edge 
& Steward 1991; David et al. 1993; Wu et al 1999; Helsdon & Ponman 2000; Xue & Wu 
2000; Croston et al. 2005). Moreover, the X-ray background in the soft band given by 
similarity is found to be much higher than the observed upper limit. In order to solve these 
problems, various models of non-gravitational heating for baryonic gas have been proposed 
(e.g. Valageas & Silk 1999; Tozzi & Norman 2001; Voit et al 2002; Zhang & Pen 2003; Xue & 
Wu 2003). The common goal of these models is to violate the similarity by heating baryonic 
gas before it falls into gravity potential wells. The amount of heating is of the order of 1 
keV nucleon" 1 (Pen 1999; Wu et al 1999). 

However, most non-gravitational heating models assume that the similarity will still 
hold without the non-gravitational heating sources. This point, as pointed out above, is 
actually inconsistent with the dynamics of the cosmic baryonic gas in the non-linear regime. 
Baryonic gas, either in high- or in low-density regions, will be heated when Burgers turbu- 
lence develops, regardless of whether non-gravitational heatingi is added. Typical shocks of 
the Burgers turbulence can convert the kinetic energy of a baryon fluid with a speed of a 
few hundred km s _1 into thermal energy (He et al. 2004; Kim et al. 2005), which is of the 
order of 1 keV nucleon -1 . Therefore, the thermodynamical effect of the Burgers turbulence 
should be substantial on the L x — T relation and soft X-ray background. We need, at least, 
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to estimate the effect of the hydrodynamic heating on the X-ray emission before considering 
extra non-gravitational heating sources. 

The paper is organized as follows. §2 describes the features of the baryonic gas as a 
Burgers fluid with samples of hydrodynamic cosmological simulation. §3 presents the analysis 
of the relationships between X-ray luminosity and mass density, temperature, and velocity 
dispersion. The soft X-ray background radiation is addressed in §4. Finally, the conclusions 
and discussion are given in §5. 



2. Baryonic gas and Burgers turbulence 

2.1. Gravitational clustering and Burgers fluid 

The baryonic gas is generally assumed to be an ideal fluid satisfying the hydrodynamic 
equations 

% + -V-(l + *)v = (1) 
ot a 

^ + (v-V)v = -iv P -V0 (2) 
ot p 

d & 1 1 1 
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where p, v, S, and p are the mass density, peculiar velocity, energy density, and pressure of 

the gas respectively. The cosmic factor a describes the cosmic expansion. The term A rac j in 

Eq.(3) is given by radiative heating and cooling of the baryonic gas per unit volume. The 

gravitational potential in Eqs.(2) and (3) is given by 

V 2 = 47rGa 2 p dm 5 dm , (4) 

where Pdm(^) and <5d m are, respectively, the mean mass density and density contrast of the 
perturbations of dark matter. Here we assume that the gravitational potential is only 
produced by the dark matter mass perturbation. 

During the process of gravitational clustering, there are two types of modes for the 
perturbations of the peculiar velocity field. For growth mode, the velocity field is curl- 
free, and the vertical mode is generally decaying with time. Therefore, one can define a 
velocity potential <p by v = — (l/a)V</?. To sketch the gravitational clustering, we consider 
the case that all the thermal processes are approximated by the polytropic relations p oc p 7 , 
T oc p 7 " 1 , or T = T^l + tf) 7-1 , where <5(x, t) = [p(x, t) —p(t)]/p(t) is the baryon mass density 
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perturbation. In this case, the equation of velocity potential, Eq.(3), can be approximately 
rewritten as 

!__L (v ^_ilvW, (5, 

where the coefficient vja 2 = 1/kj acts like a viscosity due to thermal diffusion and is char- 
acterized by the Jeans length k\ = (a 2 /t 2 )(/im p /7/c B T ). 

Eq.(5) is a stochastic-force-driven Burgers equation, which contains two scales: the 
dissipation length or the Jeans length 1/kj and the correlation length of the gravity potential 
0, r c . The intensity of can be quantified by the density contrast of dark matter 5 dm . A basic 
feature of the Burgers equation is that turbulence will develop in the fluid if the following 
condition holds (e.g. Feng et al 2003): 

(hr c )^(5l m )^ > 1, (6) 

which corresponds to the condition that the Reynolds number must be larger than 1 for the 
turbulence of incompressible fluid. Therefore, the quantity on the left-hand side of Eq.(6) 
plays the same role as the Reynolds number in the original fluid (Lassig 2000). 

Burgers turbulence is qualitatively different from the turbulence of incompressible fluid. 
The latter generally consists of vortices on various scales, while the former is a collection 
of shocks. These features arise because, for growth modes, the fluid is potential and the 
velocity field is irrotational. If v — > 0, the velocity field in the Burgers prescription acquires 
singularities due to the discontinuities caused by strong shocks. The nonlinear feature of the 
velocity field of baryonic gas can be understood as a field consisting of these shocks. 



2.2. Samples of hydrodynamic simulation 

In the numerical calculation, we still use the hydrodynamic equations (l)-(3) to model 
the baryonic gas, because the decaying modes will disappear automatically, and the result 
is just on the growth modes. Since the field of baryonic gas in the nonlinear regime consists 
of strong and weak shocks, the hydrodynamic cosmological simulation should be able to 
capture shocks and calculate the thermal properties of the baryon fluid in front of and 
behind shocks. Therefore, an optimal simulation scheme has to satisfy two conditions: (1) it 
is effective in capturing shock and discontinuity transitions, and (2) it accurately calculates 
piecewise smooth functions with a high resolution. Condition (1) is obvious in studying a 
field consisting of shocks. Condition (2) is important in calculating the shock heating, which 
depends on the difference between gaseous dynamical and thermal properties of pre- and 
postshock regions. 
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We take the WIGEON (Weno for Intergalactic medium and Galaxy Evolution and 
formatiON) code (Feng et al. 2004), which is a cosmological hydrodynamic/iV-body code 
based on the Weighted Essentially Non-Oscillatory (WENO) algorithm (Harten et al. 1986; 
Liu et al. 1994; Jiang & Shu 1996; Shu 1998; Fedkiw, Sapiro & Shu 2003; Shu 2003). The 
WENO algorithm is an Eulerian approach. It has been applied to hydrodynamic problems 
containing strong shocks, complex structures, and turbulence. It has also been applied to the 
Burgers equation (Shu 1999). The WIGEON code has passed the necessary tests including 
the Sedov blast wave and the formation of the Zel'dovich pancake (Feng et al. 2004). It 
is successful in showing the features of the Burgers fluid with the velocity field of cosmic 
baryonic gas. That is, when the Reynolds number is high, the velocity field consists of 
an ensemble of shocks with the following features: (1) the probability distribution function 
(PDF) of velocity is asymmetric between acceleration and deceleration events and (2) the 
PDF of velocity difference Av = v(x + r) — v (x) satisfies the scaling relation for a Burgers 
fluid (Kim et al. 2005; He et al. 2005). 

In order to accurately calculate the shock heating, the resolution of the simulation should 
be less than the thickness of the shock, which is of the order of the Jeans diffusion, which is 
greater than 0.2 h^ 1 Mpc for redshifts z < 2 (Bi et al. 2003). On scales less than the Jeans 
diffusion, shock-caused discontinuity is small, and the heating effect is weak. We perform 
three sets of simulations: two (samples A and B) are in a periodic cubic box with a size of 
100 hr 1 Mpc, a 512 3 grid and an equal number of dark matter particles, and another one 
(sample C) is in a cubic box with a size of 25 hr x Mpc, a 192 3 grid and an equal number 
of dark matter particles. The sizes of the grids are 100/512 = 0.20 and 25/192 = 0.13 hr x 
Mpc, respectively. Sample C is the same as that used in the analysis of the Burgers fluid 
features of cosmic baryonic gas (Pando et al 2004; Kim et al 2005). 

We use the standard ACDM model, which is specified by the matter density parameter 
Q m = 0.27, baryonic matter density parameter f2 b = 0.044, cosmological constant fl\ = 0.73, 
Hubble constant h = 0.71, and mass fluctuation erg = 0.84 within a sphere of radius 8 h" 1 
Mpc. The transfer function is calculated using CMBFAST (Seljak & Zaldarriaga 1996). For 
sample C, we use the same parameters but f2b — 0.026. We use the cloud-in-cell method for 
mass assignment and interpolation and adopt the seven-point finite difference to approximate 
the Laplacian operator. The simulations start at the redshift z = 49, and the results output 
at redshifts z=2, 1, 0.5, and 0. 

For the two simulations of the 100 h^ 1 Mpc cubic box, one considered the metal abun- 
dance (sample A), and one used the primordial composition (sample B). For sample C, we 
used the primordial composition, for which the atomic process of H and He (X=0.76,F=0.24) 
is calculated as in Theuns et al. (1998). For sample A, the metal cooling and metal line 
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emission is calculated by the phenomenological method (1) assuming a uniform evolving 
metallicity Z = 0.3Z Q (t/t ), where to is the present universe age, and (2) computing the 
cooling function using the table of Sutherland & Dopita (1993). An evolving background 
UV spectrum calculated by Haardt & Madau (1996) is used. With the comparison among 
samples A, B, and C, one can estimate the effects of the size of the simulation box and the 
metal abundance. 

A common problem of hydrodynamic simulation with Eulerian variables is that it cannot 
describe compact objects on scales less than the size of the grid. Therefore, we do not use 
these samples to study rich clusters. On the other hand, the samples with Eulerian variables 
make it easy to reach low-density regions. They are suitable for analyzing X-ray emission of 
baryonic gas from background and weakly clustered structures, such as groups. 

2.3. Shock heating and similarity breaking 

Shock heating is well known in the heating problem of groups, clusters (Cavaliere et 
al 1997), and large-scale structures (Miniati et al 2004; Kang et al. 2004). The scenario 
of Burgers turbulence shows that, in nonlinear regime, shocks or complex structures of the 
velocity field form around massive halos, as well as regions with moderate and even low mass 
density. Therefore, it heats the gas in various density fields. Figure 1 presents the relations 
between the temperature and mass density of baryonic gas at redshifts z =2, 1, 0.5, and 0.0. 
Each panel contains 19,200 data points randomly drawn from sample C. Figure 1 shows that 
the gas is multiphase. It contains the phase of tightly correlated T and p in the region of 
T < 10 4,5 K and 0.01 < pi gm < 2 and also a phase of the scattered distribution of T (> 10 5 K) 
with respect to a given p. Comparing with other numerical simulation results (e.g., Ryu et 
al. 1993; Katz, Weinberg, & Hernquist 1996; Theuns et al. 1998; Dave et al. 1999; Valageas, 
Schaeffer, Silk 2002; Springel & Hernquist 2002), the shock-heated gas is more significant 
in the moderate- and low-density (p igm ~ 1) regions. This is because the WENO code is 
effective in capturing shocks in both high- and low-density regions. 

A dynamical effect of the Burgers turbulence is the discrepancy between the fields of 
baryonic gas and dark matter. Figure 2 shows two-dimensional contours of the baryonic gas 
density p igm and dark matter density pdm, both of which are in units of their corresponding 
average density p igm and pd m - One can clearly see from Figure 2 that the density of baryonic 
gas (right) is not simply proportional to that of dark matter (left). The size of the over- 
density region p; gm > 1 is generally larger than the dark matter counterpart. Therefore, the 
clumping of baryonic gas is generally weaker than that of dark matter. 
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Figure 3 shows the same slice as in Figure 2, but gives the two-dimensional contours of 
temperature and the ratio between the densities of baryonic matter and dark matter. The 
left panel shows that the contours of Pi gm /Pdm > 1 are located outside the massive halos. 
Contrarily, the center of the halos has only a low baryon fraction Pi gm /Pdm < 1- It is more 
interesting to see from the right panel that the shape of high-temperature regions is very 
different from high-density regions of either dark matter or baryonic matter. In Figure 3, the 
high-density regions of p igm or p dm consist of filaments and knots, while the high-temperature 
regions are isolated spots. Therefore, temperature is not simply proportional to p igm or p dm , 
but multiphase. 

In summary, with respect to the similarity model, the non-linear evolution of cosmic 
baryon fluid leads to (1) a shallower distribution of baryonic gas, (2) a lower density pi gm 
inside the dark matter halo, and (3) a higher temperature outside dark matter halos. Gen- 
erally, various non-gravitational heating models try to realize the above-mentioned features. 
Figures 1-3 show, however, that all of these features are naturally yielded when baryonic gas 
is undergoing a Burgers turbulence evolution. 



3. X-ray emission of baryonic gas from groups of galaxies 

3.1. Relationship between X-ray luminosity and mass density 

For the baryonic gas in each cell of the grid with temperature T > 10 5 K, we can 
calculate the X-ray luminosity by the bremsstrahlung process assuming that all H and He 
atoms are ionized, and then obtain a field of luminosity L x . In Figure 4, we plot, respectively, 
the relationships between the X-ray luminosity and mass density of baryonic gas at redshifts 
£=0.0, 0.5, 1.0, and 2, which are given by 9931, 9119, 7314, and 4864 data points randomly 
selected from sample C. Here we consider only the bremsstrahlung radiative process, without 
considering the effect of metal line emission, which may lead to an uncertainty of no more 
than about 20%. 

Figure 4 shows that for z < 2, the L x — p igTn distribution is weakly dependent on redshift, 
and X-ray luminosity is correlated with the density p igm . Considering the adiabatic "equation 
of state" T oc p 2/3 (e.g. He et al 2004), we then obtain L x oc pl gni T l l 2 oc p 23 . However, the 
L x -Pigm distribution is actually scattered due to the multiphase property of baryonic gas. 
For a given pi gm , the scatter of baryonic gas temperature is of the order of 10 3 , i.e. from 10 5 
to 10 7 K. Therefore, the scatter of L x for a given p igm can be as large as a factor of 10 L5 . 

Table 1 gives the mean X-ray luminosity for sample C from regions with temperature 
10 6 < T < 10 7 K and dark matter density pd m < 100, 50, and 10. It shows that the X-ray 
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luminosity is mostly produced from the epoch of z < 1. At z — 1, about 25% of the X-ray 
emission is from moderate and low clustering regions pd m < 50, while at z — 0, only about 
10% is from moderate- and low-clustering areas. 

Figure 5 presents the L x — pd m relationship for the z = sample, but the fields of p dm 
and T are decomposed to the cells on scales of 0.52, 1.04, and 2.08 hr x Mpc. We use the 
scaling function of the discrete wavelet transform (DWT) to perform the decomposition. 
The DWT decomposition does not cause false correlation (Fang & Thews 1998). Therefore, 
L x of Figure 5 is the total luminosity from the cells on scales of 0.52, 1.04, and 2.08 h^ 1 
Mpc, and pd m is the mean density of dark matter in these cells. 

The cells on scales of 0.52, 1.04, and 2.08 h' 1 Mpc with density p dm > 50 have total 
masses > 0.27,2.18, and 17.4 x 10 12 M Q , respectively. These masses and sizes correspond 
to groups or clusters. It has been shown that the DWT-identified cells on scale of 1.5 
h~ x Mpc with high density are statistically the same as clusters identified by a traditional 
method, such as the friends-of- friends algorithm (Xu et al. 1998). Moreover, the cells with 
high density are found to be virialized or quasi- virialized (Xu et al 2000). Therefore, the 
DWT-identified cells with p dm > 50 on a scale of 1.04 fr 1 Mpc or p dm > 100 on a scale of 
0.52 h -1 Mpc give an ensemble of groups and clusters. This can be applied to estimate the 
statistical properties of galaxy groups. The cells with high L x simulate the behavior of the 
X-ray emission of galaxy groups. 

Generally, the mean of the luminosity L x of Figure 5 becomes larger with the increase 
of the cell scale. This is trivial, because the volume of a 1.04 (2.08) h~ l Mpc cell is 8 times 
larger than that of a 0.52 (1.04) h" 1 Mpc cell, and so on. An interesting feature shown in 
Figure 5 is that the scatter of the L x — T distributions does not reduce with the increase 
of cell scale, even when the smoothing scale of 2.08 h^ 1 Mpc is larger than the largest size 
of virialized collapsed objects. Therefore, the scatter of luminosity should be not only from 
collapse, but dependent on the multiphase property of the temperature field of baryonic gas. 



Table 1. Mean intensity of X-ray luminosity (10 6 < T < 10 7 K) 



redshift 


total 


Pdm < 100 Pdm < 50 


Pdm < 10 






(in unit of 10 43 ergs s" 1 ) 




z = 2.0 


9.10xl0~ 8 


4.12xl0~ 8 2.69xl0~ 8 


5.39xl0~ 9 


z = 1.0 


1.12xl0~ 7 


4.41xl0~ 8 2.77xl0~ 8 


5.04xl0~ 9 


z = 0.0 


2.74xl0~ 7 


5.27xl0~ 8 3.28xl0~ 8 


5.63xl0- 9 
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3.2. Relationship between X-ray luminosity and temperature 

Figure 6 presents the distribution of L x versus T for the simulation sample A at z = 
on a scale of 0.78 h~ x Mpc. The observed L x — T distribution of galaxy groups is also shown 
in Figure 6. The data are taken from Helsdon & Ponman (2000), Xue & Wu (2000), and 
Croston et al (2005). The three observed samples are partially overlapped. If a group is 
listed in two or three samples, we only use the data of the latest one. 

Figure 6 shows that the simulated results are in good agreement with the observed 
L x — T distribution of groups. First, the simulated points of high L x and T are located in 
the same areas as in the observed L x — T distribution. As mentioned in §3.1 above, the high 
L x cells are virialized or quasi- virialized objects with masses of groups. Second, the observed 
L x — T distribution is connected smoothly and continuously with simulation points with 
lower L x and T. The cells with lower L x have lower density pi gm or pdm- Therefore, Figure 
6 strongly indicates that the X-ray emission from baryonic gas related to groups underwent 
the same dynamical and thermodynamical evolution as that of the diffuse clouds, which have 
weaker X-ray emission than groups. 

Similar to the distribution of L x versus p; gm and pdm (Figures 4 and 5), the L x — T 
distribution is largely scattered. This scatter cannot be explained by the Jeans diffusion or 
another Gaussian noise process. Figure 6 shows an upper envelope in the L x —T distributions, 
which corresponds to the maximum gas density pi gm for a given temperature T. The upper 
envelope can be approximately fitted by a power law L x oc T 2 5 . Thus, we can conclude that 
if we try to fit the L x — T relationship of groups by a power law L x oc T aLT , the index a LT 
should be > 2.5. This result is consistent with the measured o^x for samples of groups, 
which are 3.6 < a < 8.2 (Helsdon & Ponman 2000), 2.1 < a < 5.7 (Xue & Wu 2000), and 
2.7 < a < 4.1 (Croston et al. 2005). 

Figure 7 is the same as Figure 6, but for sample C. The DWT variables are on scales 
of 0.52 and 1.04 hr x Mpc. The simulation data in Figure 7 are also basically consistent 
with the observed result. Therefore, the effect of metal cooling on the L x — T relation is 
small. As we know, the major effect of the metal cooling is to form a condensation branch 
in the temperature-density distribution. The temperature of the condensation branch is 
low, and therefore, these are not strong X-ray sources. It is interesting to see that Figure 7 
shows a lower envelope with a power law L x oc T~ 5 . This power law is about the same as 
L x oc T~ 4 8 given by the hydrodynamic simulation of Dave et al (2002). However, our result 
of the L x — T distribution is scattered, i.e. fills in the range between the lower and upper 
envelopes. Observed samples also show the scatter and are consistent with our simulation 
sample. 
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Figure 8 presents the L x — T distribution for sample C, but decomposed on physical 
scales of 0.69 hr x Mpc at redshift z = 0.5 and 0.5 hr x Mpc at z — 1. The observed points 
are the same as those in Figures 6 and 7. We see that the simulated L x — T distribution is 
almost z independent in the range z < 1. This is consistent with a weak evolution scenario 
of temperature and entropy fields in the range z < 2 (He et al. 2004). Thus, we may predict 
that the L x — T distribution of groups at z ~ 1 should not be very different from z = 0. 

3.3. Relationship between X-ray luminosity and velocity dispersion 

The velocity dispersion in a cell can be estimated by the peculiar velocity difference 
within this cell. We use the DWT variables to calculate the velocity dispersion. The details 
of this algorithm are given by Yang et al (2001) and Kim et al. (2005). This method is 
effective in estimating the velocity dispersion of both virialized systems and quasi- virialized 
systems (Xu et al. 2000). 

Figure 9 plots the distribution of X-ray luminosity L x versus velocity dispersion of dark 
matter a^m for simulation sample C decomposed in cells on a scale of 1.04 h^ 1 Mpc. The 
observed L x - a points of Mulchaey & Zabludoff (1998) are also shown in Figure 9. Many 
objects in this sample are actually rich clusters, which are outside of the simulated data 
points. Nevertheless, the observed points with low L x and a are clearly support to the 
simulated result. 

The distribution of L x versus <7d m in Figure 9 is also substantially scattered. However, 
it does dot have a very clear upper envelopes but one can still see a less clear upper envelop, 
which can be roughly fitted by a power law L x oc a 21 . That is, if we try to describe the 
scattered L x — a distribution by a power law L x oc a aLV , the index aiv should be larger 
than 2.1. This is consistent with observational results aiv > 2.3 (Helsdon & Ponman 2000; 
Xue & Wu 2000). The large scatter of L x versus a is caused by the large scatter of velocity 
versus density of both baryonic matter and dark matter. 

4. Soft X-ray background 

4.1. Mean and PDF of soft X-ray background 

As the hot gas with high temperature (T > 10 7 K) is always identified as clusters, we 
calculate the X-ray background intensity only from baryonic gas with temperature 10 5 — 10 7 K, 
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which is given by an integral 

1 c f Zr dz 
F ^ ) = T,Wj e ([l + ^.,,) (1 + 2)1[SJm(1 + 2)3 + sjA]1/2 , (7) 

where F(vq) is in units of ergs s _1 cm -2 Hz -1 sr _1 and e(uo,z) is the volume emissivity of 
X-ray photons with local frequency u at redshift z. The emissivity e(z/ , z) is calculated 
using Raymond-Smith model (Raymond & Smith 1977). In the calculation of the X-ray 
background intensity, both the samples with and without metal cooling are considered. 

The soft X-ray background mostly comes from hot clouds at z < 2. It is sufficient to 
produce maps by taking the integral of Eq.(7) from z = to 6. The algorithm of mapping 
is as follows. The simulation data are output every time when light crosses the simulation 
box with a size of 100 h _1 Mpc. For each output, we take the integral in Eq.(7) along one 
random chosen axis and produce a map of X-ray flux F at the redshift z of that box. After 
the simulation reaches z = 0, we have 60 maps of F at redshifts corresponding to a comoving 
distance at n x 100 h^ 1 Mpc, where n = 0, ...59. One map of the background can then be 
obtained by (1) arranging a set of 60 two-dimensional maps with a randomized center, (2) 
taking a constant angular projection for each map, and (3) superposing the contributions of 
the 60 maps. 

A typical two-dimensional map of the soft X-ray background from regions with 10 5 < 
T < 10 7 K based on sample A is shown in Figure 10. The angular size of the map is 1° x 1°, 
and its angular resolution is 0'.06. Figure 11 gives the corresponding PDF of the intensity of 
the X-ray background in the band 0.5 — 2 keV. For sample A, the mean X-ray background 
intensity is (1.87±0.63) x 10~ 12 ergs s _1 cm" 2 deg -2 , where the 1 a error is from the variance 
of 50 two-dimensional maps. For sample B, we have (1.73±0.63) x 10~ 12 ergs s -1 . Therefore, 
the effect of metal cooling has to be less than 10%. 

Table 2. Mean intensity of soft-band X-ray background 
p dm Mean intensity (ergs s" 1 cm -2 deg -2 ) 



^10 5 < T < 10 7 K) Primordial abundance Metal abundance 

total 1.73±0.63 x 10~ 12 1.87±0.63 x 10~ 12 

p dm < 100 2.90±0.62 x 10- 13 2.86±0.86 x 10~ 13 

Pdm < 50 1.27±0.29 x 10~ 13 1.23±0.38 x 10" 13 

p dm < 10 1.55±0.41 x 10" 14 1.52±0.63 x 10~ 14 
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4.2. Density dependence of soft-band X-ray background 

It is usually believed that most X-ray background is attributed to resolvable discrete 
X-ray sources, such as active galactic nuclei (AGNs) and truly diffuse extragalactic baryonic 
clouds. As shown in §3, the dynamical and thermodynamical properties of baryonic clouds 
are continuously varying from low-density (diffuse) regions to high-density regions (collapsed 
and virialized objects). There is no clear criterion to distinguish the X-ray emission of truly 
diffuse baryonic gas from that of resolvable sources. We study this problem with the density 
dependence of the soft X-ray background. 

We integrate Eq.(7), but only calculate X-ray emission from baryonic clouds with tem- 
peratures 10 5 < T < 10 7 K and in regions with p dm <100, 50, and 10. The results for sample 
A are shown in Figure 12, which contains three 1° x 1° maps corresponding to pdm <100, 
50, and 10. The mean intensities of the soft-band X-ray background are listed in Table 2. 
We can see that less than 17% (7%) of the total soft X-ray background of 10 5 < T < 10 7 
comes from clouds located in regions of dark matter with density contrast pd m < 100 (50). 
As mentioned in §3.2, the cells having pd m > 100 or > 50 generally contain collapsed and 
virialized structures. From Figure 2 we can also see that the structures with pd m > 100 or 50 
have only a very small volume fraction. These regions are discretely distributed. Therefore, 
they should be resolvable. Thus, one can conclude that the contribution of truly diffuse 
baryonic gas to the soft-band X-ray background is no more than about 2.90 x 10~ 13 ergs s -1 
cm -2 deg -2 . This result can also be seen from Table 1. Using data of the XMM — Newton 
Lockman Hole observation, Worsley et al (2005) show that in the soft-band (<2 keV), more 
than 90% of the X-ray background can be resolved. This result is consistent with Table 
2. The simulation done by Croft et al (2001) found the X-ray background of the warm-hot 
intergalactic medium (WHIM) to be 4.15 x 10 13 ergs s -1 cm -2 deg -2 . This is a little larger 
than our result. This difference is probably attributable to the powerful shock-capturing 
ability of the WENO code. 

It is also interesting to note from Table 2 that for the case of the total region, the 
sample considering the metal abundance (sample A) has a slightly higher mean intensity of 
the soft X-ray background than that without metal cooling (sample B), while for all other 
cases (pdm < 100) sample A gives a slightly lower mean intensity than that of the primordial 
abundance. This is because metal cooling can weaken the turbulence heating. The gas with 
metal cooling easily falls into gravity wells (p dm > 100). On the other hand, in the case of no 
metal cooling more gas will remain in less clustered regions and give higher X-ray emission 
from clouds with pd m < 100. In all cases, the difference between samples A and B is small. 
Thus, the thermodynamical effect of metal cooling is much less the turbulence heating. 
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4.3. Temperature dependence of soft X-ray background 

Similarly, we study the temperature dependence of soft X-ray background. We do the 
integral in Eq.(7), but only calculate X-ray emissions from baryonic clouds with temperatures 
less than 10 6 ' 5 K. A 1° x 1° map of the soft X-ray background from clouds of T < 10 6 ' 5 of 
sample A is shown in Figure 13, and the corresponding PDF of the intensity of the X-ray 
background is given in Figure 14. The mean intensity is (9.38 ± 3.06) x 1CT 14 ergs s" 1 cm" 2 
deg~ 2 for sample A and (8.0 ±2.03) x 10~ 14 ergs s _1 cm -2 deg -2 for sample B. Therefore, no 
more than 5% of the total soft X-ray background of 10 5 < T < 10 7 K comes from clouds with 
temperature T < 10 6 ' 5 K. Most of soft X-ray background is attributed to baryonic clouds 
with high temperature (T > 10 6 5 ). 

5. Discussion and conclusions 

The evolution of baryonic gas can be roughly divided into three stages. The first is the 
linear stage. It can be simply described by the similarity between the baryonic and dark 
matter. When the non-linear evolution takes place, Burgers turbulence develops. Baryonic 
gas will decouple with dark matter, and becomes multiphase. The kinetic energy of baryon 
fluid will be dissipated by the shocks of turbulence. Some cosmic baryonic gas finally falls 
into massive halos and evolves in the thermal equilibrium state in gravity wells. Some 
heated baryonic gas remains outside of gravity wells. The Burgers turbulence heating is not 
uniform. A uniform heating would contradict with the Lya forest, which is produced by 
hydrogen clouds with a temperature of about 10 4 -10 5 K. Therefore, baryonic gas undergoing 
the Burgers turbulence evolution must be multiphase. In the stage of the Burgers turbulence, 
the dominant component in volume fraction is the gas in the low-temperature phase (He et 
al. 2004). 

We show that the X-ray emission of baryonic gas can be well modeled by the above- 
mentioned scenario. That is, the observed distributions of L x versus T and L x versus a 
of galaxy groups can be reproduced by the simulation sample of the ACDM model without 
adding extra heating sources. We also find that almost all of the soft-band X-ray background 
radiation is from clustered regions, and the contribution of diffuse gas with p^ m < 50 is 
negligible. This is also consistent with current X-ray observations. These results show that 
the Burgers turbulence heating alone seems to be enough to account for the basic features 
of X-ray emission of baryonic gas in the universe. 

We should also point out that star formation and its feedback on the baryonic gas are not 
considered in our simulation. Roughly, there are two types of feedback: (1) photoionization 
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heating by the UV emission of stars and AGNs and (2) injection of hot gas and energy by 
supernova explosions, or other sources of cosmic rays. Actually, the photoionization heating 
can be properly considered, if the UV background is adjusted by fitting the simulation with 
the observed mean flux decrement of QSOs' Lya absorption spectrum (Feng et al. 2003). The 
significant effect of injecting hot gas and energy by supernovae is mostly on dwarf galaxies. 
Recent Lya observations of protoclusters (Adelberger et al. 2003) suggest that AGN heating 
does not drastically affect the gas in clusters. Therefore, this heating mechanism may not 
be strong enough to change the basic features of the multiphase scenario given in this paper. 
Moreover, the shocks may lead to the electron temperature T e being lower than the ion 
temperature T ; . Since the mean of the ratio T-JT C is generally less than 3 (Fox & Loeb 1997; 
Takizawa 1998; Yoshida et al. 2005), this effect may lead to an uncertainty of the luminosity 
L x that is less than a factor of 2. Therefore, all conclusions should basically be held. 
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Fig. 1. — Temperature T vs. baryon density p igm for sample C at redshifts 2, 1, 0.5, and 0, 
where p igm is in units of the mean baryonic matter density p igm . Each panel contains 19,200 
data points randomly selected from sample C. 
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Fig. 2. — Density contour plots of dark matter (left) and baryonic gas (right) for a slice of 
0.26 h^ 1 Mpc thickness at z — 0. The solid contours encompass the over-dense regions with 
p = e*/ 2 , i — 0,1, 2... (p is normalized to 1), while the dotted lines represent the boundaries 
of the under-dense regions with p = e~ % l 2 , % = 1,2.... 
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Fig. 3. — Contour plots of baryon to dark ratio Pi gm /pdm {left) and temperature T {right) 
for the same slice as in Figure 2. The solid contours represent, respectively, the regions with 
Pigm/Pdm = e i/3 ,i = 0,1,2,..., and T = e i/2 x 10 5 K, i = 0,1,2,..., while the dotted lines 
represent Pi gm /Pdm = e~ 1 / 3 , i = 1,2, and T = e -1 / 2 x 10 5 K, % — 1, 2, regions. 
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Fig. 4. — X-ray luminosity L x vs. baryon density p; gm for sample C, where p; gm is in units 
of the mean baryonic matter density Pi gm . The panels at z—0.0, 0.5, 1.0, and 2 show, 
respectively, 9931, 9119, 7314, and 4864 data points randomly selected from sample C in 
each redshift, and with temperature T > 10 5 K. 
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Fig. 5. — X-ray luminosity L x vs. dark matter density pdm for sample C at redshift z = 
and decomposed on scales of 0.52, 1.04, and 2.08 h~ 1 Mpc. Here pdm is in units of the mean 
dark matter density pdm- 
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Fig. 6. — X-ray luminosity L x vs. temperature T for samples A at redshift z = with 
the DWT decomposition on scales of 0.78 h~ x Mpc. The observed data are adapted from 
Helsdson & Ponman (2000; squares), Xue & Wu (2000; circles) and Croston et al (2005; 
triangles) . 
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Fig. 7. — Same as Fig.6, but for sample C at redshift z = with the DWT decomposition 
on scales of 0.52 (top) and 1.04 h" 1 Mpc (bottom). 
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Fig. 8. — Same as Fig. 7, but for redshift z = 0.5 (top) and z = 1 (bottom). The physical 
scale of the DWT decomposition is 0.69 h^ 1 Mpc. 
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Fig. 9. — X-ray luminosity L x vs. velocity dispersion of dark matter a^m for sample C at 
redshift z = with the DWT decomposition on scales of 1.04 h~ x Mpc. The observed data 
are adapted from Mulchaey & Zabludoff(1998). 
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Fig. 11. — Probability distribution function of X-ray intensity in the 0.5-2 keV band for the 
same samples as in Fig. 10. 
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Fig. 12. — Same as Fig. 10, but for baryonic clouds located in regions with dark matter 
density pdm < 100 (top), 50 (middle) and 10 (bottom), where pdm is in units of the mean 
dark matter density pdm- 
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Fig. 14. — Same as Fig. 11, but for the x-ray emission from clouds with temperature in range 
10 5 < T < 10 6 5 K (the same samples as in Fig. 13). 



